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Abstract 

The dynamics of the second order rational difference equation Zn+i = . ^ T. ^ 

A + BZn + CZn-l 

with complex parameters and arbitrary complex initial conditions is investigated. 

In the complex set up, the local asymptotic stability and boundedness are studied 
vividly for this difference equation. Several interesting characteristics of the solu¬ 
tions of this equation, using computations, which does not arise when we consider 
the same equation with positive real parameters and initial conditions are shown. 

The chaotic solutions of the difference equation is absolutely new feature in the 
complex set up which is also shown in this article. Some of the interesting observa¬ 
tions led us to pose some open interesting problems regarding chaotic and higher 
order periodic solutions and global asymptotic convergence of this equation. 
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1 Introduction and Preliminaries 

A rational difference equation is a nonlinear difference equation of the form 


^n+l 


“ 1 “ ^—n 

^ + Si =0 BiZn-i 


where the initial conditions Zq, z^i,..., z^k are such that the denominator never vanishes 
for any n. 

Consider the equation 


a + /^Zn + 7Zn-l 

where all the parameters and the initial conditions z^i and Zq are arbitrary complex 
number. 

This second order rational difference equation Eq. Q is studied when the parameters are 
real numbers and initial conditions are non-negative real numbers in |2]. In this present 
article it is an attempt to understand the same in the complex plane. 


Here, a very brief review of some well known results which will be useful in order to 
apprehend the behavior of solutions of the difference equation ([^. 

Let / : —)■ D where D C C be a continuously differentiable function. Then for any 

pair of initial conditions Zq, z^i G D, the difference equation 

^n+l f {,^n^ ^n—l) (2) 


with initial conditions Z-i,Zo G D. 

Then for any initial value, the difference equation (1) will have a unique solution {zn}n- 
A point ^ G D is called equilibrium point of Eq.(|^ if 

f{z,z) = z. 

The linearized equation of Eq. (|^ about the equilibrium ^ is the linear difference equation 


where for i = 0 and 1. 


^n+l T (^l^n—1 y ^ 0, 1, . . . 

df , 


The characteristic equation of Eq.(2) is the equation 


( 3 ) 
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— oqA — Oil — 0. (4) 

The following are the briehngs of the linearized stability criterions which are useful in 
determining the local stability character of the equilibrium z of Eq. (@, n. 

Let z be an equilibrium of the difference equation Zn+i = fizn, Zn-i)- 

• The equilibrium z of Eq. (2) is called locally stable if for every e > 0, there 
exists a 5 > 0 such that for every zq and z_i G C with \zo — z\ + |z_i — z| < <5 we 
have \zn — z\ < e for all n > — 1 . 

• The equilibrium z of Eq. (2) is called locally stable if it is locally stable and if 
there exist a 7 > 0 such that for every zq and z_i G C with I^q — z\ + |^_i — z| < 7 
we have lim^^oo Zn = z. 

• The equilibrium 7 of Eq. (2) is called global attractor if for every zq and z_i G 
C, we have lim^^oo Zn = z. 

• The equilibrium of equation Eq. (2) is called globally asymptotically stable/fit 
is stable and is a global attractor. 

• The equilibrium 7 of Eq. (2) is called unstable if it is not stable. 

• The equilibrium 7 of Eq. (2) is called source or repeller if there exists r > 0 
such that for every zq and Z-i G C with \zq — z\ + \z-i — z\ < r we have \zn — z\ > r. 
Clearly a source is an unstable equilibrium. 

Result 1.1: (Clark’s Theorem) The sufficient condition for the asymptotic stability 
of the difference equation ( 1 ) is 

|®o| + |®i| < 1 


2 Difference Equation and Its Transformed Forms 


The following difference equation is considered to be studied here. 

a + /3Zn + 7^n-l „ , 

Zn+1 = . , p -, n = 0, 1, . . . 

A + BZn + CZn-1 

where all the parameters are complex number and the initial conditions z^i and Zq are 
arbitrary complex numbers. 


We will consider three different cases of the Eq.(|^ which are as follows: 
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a 


2.1 The case P = j = 0 : Zn+i 


A+BZn + CZn-l 


By the change of variables, Zn = ^Wn, the difference eqnation Zn+i 
rednced to the difference eqnation 


_ a _ 

A-\-BZn-\-CZn — l 


^n +1 


where p = ^ and q = 


, n = 0, 1, . . . 

1 + pWn + qWn-1 


( 6 ) 


2.2 The case a — j — 0 : Zn+i — A+Bzn+Czn-i 

By the change of variables, Zn = pWn, the difference eqnation Zn+i 
reduced to the difference equation 


__ 

A-\-BZn-\-C Zn — \ 


^n +1 


where P = ^ ^ ~ 


Wr, 


1 +pWn + qWn-l 


,n = 0,1, 


(7) 


2.3 The case a = 13 = 0 : z„+-i = ^+b':Zcz„., 

By the change of variables, Zn = pWn, the difference equation Zn+i 
reduced to the difference equation 


Wn+l = 


'^n—1 


p + qWn + Wn-l 


,n = 0,1, 


where p = - and q = 


c 


'yZn-1 

A-\-BZn-\-CZn, — l 


( 8 ) 


Without any loss of generality, we shall now onward focus only on the three difference 
equations (6), (7) and (8). 


3 Local Asymptotic Stability of the Equilibriums 

In this section we establish the local stability character of the equilibria of Eq.Q in 
three difference cases as stated in the section 2. 

3.1 Local Asymptotic Stability of Wn+i = tt— r - 

L-\-pWn~rQ'^n—l 

The equilibrium points of Eq.(6) are the solutions of the quadratic equation 

1 

W = - -z-r 

1 + pw + qw 
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Eq.(6) has the two equilibria points Wi ^2 = ~^~ 2 (pX*q)^ respectively. 

The linearized equation of the rational difference equation(6) with respect to the equi¬ 
librium point wi = is 


4p 


-,Wn + 


4g 


Zn+l + -- ,_, ^ . , ,_. 

(-1 ^/l + Ap + Aq) (-1 -h ^/l + Ap + Aq) 

with associated characteristic equation 

4p 4g 


;Wn-i = 0,n = 0,1,... (9) 


A" + 


A + 


(-1 -h a/1 + 4p + dg) (-1 + y/l + Ap + Aq) 


= 0 . 


( 10 ) 


The following result gives the local asymptotic stability of the equilibrium wi of the Eq. 

( 6 ). 


Theorem 3.1. The equilibriums Wi = ^ Xl+q)^ of Eq. (6) is 

locally asymptotically stable if 


Ap 

+ 

Aq 

(-1 + +Ap + 4g)“ 


(-1 + ^/l + Ap + 4g)“ 


Proof. The zeros of the characteristic equation (10) has two zeros which are 
- , ^ and- , ^ . Therefore by Clark’s 

p+i/p2-4p(j+2(j(-l-2(j-|-v^l-|-4p+4ij) -p+/p2_4pq_|_2q(_i_2q+yi+4p+4(/j 

theorem, the equilibrium wi = is locally asymptotically stable if the sum of 

the modulus of two coefficients is less than 1. Therefore the condition of the polynomial 


(10) reduces to 


Ap 


(-l+Vl+4p+4gy 


+ 


4q 


( - 1 +VT+4p+4g) 


< 1 . 


□ 


The linearized equation of the rational difference equation (6) with respect to the equi¬ 
librium point W 2 = — is 


Wn+l + 


4p 


jWn + 


4g 


(l -h a/ 1 + 4p -h 4g) (l -h a/ 1 + 4p Aq) 

with associated characteristic equation 

4p Aq 


Wn-I = 0,n = 0,1,... 


_l_ y:_y _j_ 

[l + y/1 + Ap + AqY [l + y/1 + Ap + Aq)‘ 


= 0 . 


( 11 ) 

( 12 ) 


5 



Theorem 3.2. The equilibriums W 2 = of Eq. (6) is 

locally asymptotically stable if 


4p 

+ 

Aq 

{1 + ^/l + 4:p + Aqf 


(l + + Ap + Aqf 


< 1 


Proof. Proof the theorem follows from Clark’s theorem of local asymptotic stabil¬ 
ity of the equilibriums. The condition for the local asymptotic stability reduces to 


4p 


(l + V^T+ip+Iq)^ 


+ 


4q 


(l+Vl+4p+4g)" 


< 1 . 


□ 


Here is an example case for the local asymptotic stability of the equilibriums. 


For p = \ and | the equilibriums are —1.6838 -|- 1.13355i and 0.683802 — 0.1335521 
For the equilibrium 1.6838 -|- 1.13355i, the coefficients of the characteristic polynomial 
(10) are 0.775125 -|- 1.90868i and —1.90868 — 0.775125i with same modulus 2.06006. 
Therefore the condition as stated in the Theorem 3.1 does not hold. Therefore the 
equilibrium 1.6838 -|- 1.13355i is unstable. 


For the equilibrium 0.683802 — 0.133552z, the coefficients of the characteristic polyno¬ 
mial (12) are —0.224875 -|- 0.091323z and —0.091323 — 0.224875z with same modulus 
0.242711. Therefore the condition as stated in the Theorem 3.2 is hold good. Therefore 
the equilibrium 0.683802 — 0.133552i is locally asymptotically stable. 


It is seen that in the case of real positive parameters and initials values, the positive 
equilibrium of the difference equation ( 6 ) is globally asymptotically stable [2]. But the 
result is not holding well in the complex set up. 


3.2 Local Asymptotic Stability of Wn+i = tt—- 

i H-p “r ^ n — 1 

The equilibrium points of Eq.(7) are the solutions of the quadratic equation 

w 

W = - -z-r 

1 -|- pw -|- qw 

The Eq.(7) has only the zero equilibrium. The linearized equation of the rational differ¬ 
ence equation(7) with respect to the zero equilibrium is 

Wn+i = Wn,n = 0,1, ■ ■ ■ (13) 

with associated characteristic equation 

- A = 0. (14) 
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The following result gives the local asymptotic stability of the zero equilibrium of the 
Eq. (7). 

Theorem 3.3. The zero equilibriums of the Eq. (7) is non-hyperbolic. 

Proof. The zeros of the characteristic equation (14) has two zeros which are 0 and 1. 
Therefore by dehnition, the zero equilibrium is non-hyperbolic as the modulus of one 
zero is 1 . □ 

It is nice to note that in the case of real positive parameters and initials values, the 
zero equilibrium of the difference equation (7) is globally asymptotically stable for the 
parameter p > 1 [2] . But in the case of complex, the zero equilibrium is non-hyperbolic 
as we have seen the previous theorem. 


3.3 Local Asymptotic Stability of Wn+i = , ^ - 

^ P+qWn+Wn-l 

The equilibrium points of Eq.( 8 ) are the solutions of the quadratic equation 


w 

w = -z-V 

p + qw + w 

The Eq.( 8 ) has two equilibriums which are 0 and The linearized equation of the 
rational difference equation( 8 ) with respect to the zero equilibrium is 

1 , , 
Wn+l = = 0, 1, . . . (15) 

P 

with associated characteristic equation 

- - = 0 . (16) 
p 

The following result gives the local asymptotic stability of the zero equilibrium of the 
Eq. ( 8 ). 

Theorem 3.4. The zero equilibriums of the Eq. (8) is locally asymptotically stable if 
\p\ > 1 and repeller if |p| < 1. 


Proof. The zeros of the characteristic equation (15) has two zeros which are 0 and 
Therefore by dehnition, the zero equilibrium is locally asymptotically stable if 
and unstable (repeller) if - >1. Hence the required is followed. 


< 


p' 

1 

□ 


The linearized equation of the rational difference equation( 8 ) with respect to the equi¬ 
librium is 
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1+ pq p 

Wn+l = ^- Wn + 


1 + g ' 1 + g 

with associated characteristic equation 


Wn-I = 0,n = 0,1,... 


A' 


I+pq 


X 


p 


1+q l+q 


= 0 . 


(17) 


(18) 


The following result gives the local asymptotic stability of the equilibrium of the 
Eq. (8). 

Theorem 3.5. The zero equilibriums of the Eq. (8) is locally asymptotically stable if 

|1 +pq\ + Ip - 1| < |1 + g| 

Proof. The equilibrium ^ of the characteristic equation (18) would be locally asymp¬ 
totically stable if the sum of the modulus of the coefficients of the characteristic equation 


(18) is less than 1. That is by Clark’s theorem, 


i+p<? 

1+9 


+ 


P-1 

1+9 


< 1, that is 


|1 +pq\ + Ip - 1| < |1 + g| 

Here is an example case for the local asymptotic stability of the equilibriums. 


□ 


For p = 1 + I (IpI = 1.1180 > 1) and q = + * (|q'| = 1.005 > 1) the equilibriums are 0 

and —0.226244—0.248869h For the equilibrium —0.226244—0.248869*, the coefficients of 
the characteristic polynomial (18) are 0.533622 + 0.545503* and —0.203409 — 0.0330197* 
with modulus 0.763103 and 0.206072 respectively. Therefore the condition as stated 
in the Theorem 3.5 hold good. Therefore the equilibrium 1.6838 + 1.13355* is locally 
asymptotically stable. 


In the case of real positive parameters and intimal values of the difference equation (8), 
the positive equilibrium is locally asymptotically stable if p < 1 and g < 1 but in the 
complex set, it is encountered through the example above is that the equilibrium is 
locally asymptotically stable even though |p| > 1 and |g| > 1 [2]. 


4 Boundedness 

In this section we would like to explore the boundedness of the solutions of the three 
difference equations (6), (7) and (8). 

Now we would like to try to find open ball 5(0, e) G C such that if Wn G 5(0, e) and 
Wn-i G 5(0, e) then Wn+i ^ 5(0, e) for all n > 0. In other words, if the initial values wq 
and W-i belong to 5(0, e) then the solution generated by the difference equations would 
essentially be within the open ball 5(0, e). 












Theorem 4.1. For the difference equation (6), for every e > 0, ifwn and Wn-i € -B(0, e) 
then Wn+i G -8(0, e) provided 

bl > 1 + \q\ 


Proof. Let be a solution of the equation Eq.( 6 ). Let e > 0 be any arbitrary real 

number. Consider Wn, Wn-i ^ -8(0, e). We need to find out an e such that Wn+i € -8(0, e) 
for all n. It is follows from the Eq.( 6 ) that for any e > 0, using Triangle inequality for 


|Wn+l| 


1 

1 + pWn + qWn-1 


< 


1 


pWn + qWn-1 


1 

“ (IpI - kl)e 


In order to ensure that |tnn+i| < (Assuming e' = ■^) it is needed to be 

1 

(IpI - kl) ^ ^ 


That is IpI > 1 + |g|. Therefore the required is followed. 


□ 


Theorem 4.2. For the difference equation (1), for every e > 0, ifwn and Wn-i 
then Wn+i € -8(0, e) provided 


\p\ > |g| + 


1 

e 


G -8(0, e) 


Proof. Let {wn} be a solution of the equation Eq.(7). Let e > 0 be any arbitrary real 
number. Consider Wni Wn-i G .8(0, e). We need to find out an e such that Wn+i G .8(0, e) 
for all n. It is follows from the Eq.(7), that for any e > 0, using Triangle inequality for 


|Wn+l| = 


Wr, 


1 +pWn + qWn-l 


< 


i\p\ - kl)e 


Therefore, |tc„+i| < We need lyyiry to be less than e. Therefore, 


bl-kl 


\p\ ^ 1^1 + “ 


is followed. 


□ 


Theorem 4.3. For the difference equation (8), for every e > 0, ifwn and Wn-i G .8(0, e) 
then Wn+i G -8(0, e) and \p\ < 1 provided 


e < 


\P\ - 1 
|g| +1 
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Proof. Let be a solution of the equation Eq.(8). Let e > 0 be any arbitrary real 

number and |p| > 1. Consider Wn,Wn-i € -8(0, e). We need to find out an e such that 
Wn+i € -8(0, e) for all n. It is follows from the Eq.(8), that for any e > 0, using Triangle 
inequality for 

e 

\p\ - k|e-e 


kn+i| = 


Wn-l 


P + qWn + Wn-l 


Therefore, 


|Wn+l| < 


e 

\p\ - k|e-e 


In order to ensure that |tnn+i| < e, it is needed to be 


1 

\p\ - k|e-e 


< 1 


That is 


Therefore the required is followed. 


e < 


\P\ - 1 
|g| +1 


□ 


5 Periodic of Solutions 

A solution {wn}n of a difference equation is said to be globally periodic of period t if 
Wn+t = Wn for any given initial conditions, solution {wn}n is said to be periodic with 
prime period p if p is the smallest positive integer having this property. 

We shall first look for the prime period two solutions of the three difference equations 
( 6 ), (7) and ( 8 ) and their corresponding local stability analysis. 

5.1 Prime Period Two Solutions of Eq. (6) 

Let ..., 0, 0, 0, 0 7 ^ 0 be a prime period two solution of the differ¬ 
ence equation Wn+i = 1 -• Then 0 = , , (, , and tb = ,, , h —r. This 

two equations lead to the set of solutions (prime period two) except the equilib- 

f, 0.b(p-q+Jp'^+p{-2-Aq)q+q^+lqA 0.5p-0.5q-0.bJ p'^+p(-2-iq)q+q^+Aq^ 1 , 

riums as |0 ^ ^- ^,0 ^ 1 and 

_Q 0.5\/p‘^+p(-2-4q)q+q‘^+4,q3 

-; ^ ,0 ^ ^ - 

-0.5p+0.5q-0.5y'p'^+pi-2-Aq)q+q^+Aq3 ^ 9 

Let ..., 0, 0, 0, 0,..., 0 0 0 be a prime period two solution of the equation ( 6 ). We 
set 

'^n Wn—1 
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Then the equivalent form of the difference equation (6) is 

"Un+l ~ Vji 

1 

"^n+l — Z j j 

l+pVn + qUn 

Let T be the map on C x C to itself dehned by 


V 

1 

l-\-pv-\-qu 


Then 


^ ^ is a hxed point of T^, the second iterate of T. 


2 / u 


l-\-pv-\-qu 




g{u,v) 
h{u, v) 


where g{u,v) = and h{u,v) = ^ ^ p -j-^. Clearly the two cycle is locally 

asymptotically stable when the eigenvalues of the Jacobian matrix Jt 2 , evaluated at 
lie inside the unit disk. 


We have, 




where g(0,^) = and g(0,^) = 




((l+q0)(l+gV’)+p(l+V’+QV’^))^ 


and g(.#..V<) = 

l + qcP+P'iP ^ 


Now, set 




(l+q<f)+pip )2 


(l + # + pV.P (l + ,V.+ 5^) 




r \r 7 r / r 

ou ov 




((1 + g0)(l + qi^) +p{l + 'ip + q'ip‘^)y 
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In particular for the prime period 2 solution, 


-)■ 


0.b{p-q+^Jp^+p{-2-iq)q+q^+iq^'^ ^ 0.^p-0.}iq-Q.^^p2+V{-2-iq)q+q'^+Aq^ 

q{-p+q) ’ q{-p+q) 


we shall see 


the local asymptotic stability for some example cases of parameters p and q. The gen¬ 
eral form of X and A would be very complected. Consider the prime period two solution of 
the difference equation (6), 0 —)■ —0.0843748-1-0.062214500 —>■ —0.0822456 —0.0594374f 
corresponding two the parameters p —)■ 100 -|- i and g —)■ 6 -|- O.li 


In this case, |x| = 0.0735211 and |A| = 0.004075. Therefore, by the Linear Stability the¬ 
orem (Ixl < 1-1- |A| < 2) the prime period 2 solution 0 —)■ —0.0843748-1-0.062214500 —)■ 
—0.0822456 — 0.0594374f is locally asymptotically stable. 


5.2 Prime Period Two Solutions of Eq. (7) 

Let ...,0, 0, 0, 0,---, 0 0 0 be a prime period two solution of the difference 

equation Wn+i = ti —^ -• Then 0 = , , t, , and 0 = , , f, —r. This two 

equations lead to the set of solutions (prime period two) except the equilibriums as 


-)■ 


_ 1 _^ ^ _1 + VZzJ 

0.5p-0.5q-0.5^p^-q^ ’ ^ 1 


-)■ 


_ 1 _^ ^ _1 _ 

O.5p-O.5q+O.50pTr^’ ^ q pq-q^ 


Let ..., 0, 0, 0, 0,..., 0 0 0 be a prime period two solution of the equation (7). We 
set 

'^n—1 
V„. = W„, 


Then the equivalent form of the difference equation (7) is 




Un+1 — Vn 
Vn 


I +pVn + qUr 
Let T be the map on C x C to itself dehned by 


T 


u 


l-\-pv-\-qu 


Then 


0 


is a hxed point of T^, the second iterate of T. 


T2 


u 


( 


l+pv-\-qu 


\ 


l+pn + gtt 

\ i+p +9^ / 
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Ti 


g{u,v) 
h{u, v) 


where g{u,v) = 


and h{u,v) = 


l-\-pv-\-qu 


l+pv-\-qu 

able 

, , lie inside the unit disk. 

We have, 

fMi-) 

where g(,#.,i/>) = 


Clearly the two cycle is locally 


i+P i+p"+gu +‘3^' 

asymptotically stable when the eigenvalues of the Jacobian matrix Jt^, evaluated at 


Jx'^ 






!+(}((/)— 


Now, set 


(1+ 1 /’+pV’+9 (</>+'!/>++pV’2)) 


qtjj 


(i+y+pp+i}(</)+-!/'+(j</)'i/>+py2)) 


2^^2 


1 + g (0 — pip'^ 


X = ^( 0 ,^) + ^(0,^) = - r. , "" /N2 + n / / /A ! A! /2^^2 

du dv (l+g0+P0)2 (1 ^ -I- _|_ g (^0 -I- ^ -p 




A = 


(1 + g 0 + pip) {1 + Ip + pip + q{4> + Ip + q(pip + pip‘^))‘^ 
In particular for the prime period 2 solution. 


I 1 1 .» /p2_ I 

< 0 —)■- / , Ip —)■—- + - 2 “ c, we shall see the local asymptotic stabil- 

0.5p-0.5(?-0.5A/p2-g2’ ^ 9 pg-g" j’ ^ 

ity for some example cases of parameters p and q. The general form of y and A would 
very complected. Consider the prime period two solution of the difference equation (7), 
(p —)■ 0.365026 + O.263198i,0 —)■ —0.412345 + 0.131124i corresponding two the parame¬ 
ters p —)■ I J- 3i and g —| J- 57 


In this case, |y| = 0.0287948 and |A| = 0.0000431717. Therefore, by the Linear Stability 
theorem (|y| < 1 J- |A| < 2) the prime period 2 solution (p —)■ 0.365026 J- 0.26319800 —)■ 
—0.412345 J- 0.131124i is locally asymptotically stable. 


5.3 Prime Period Two Solutions of Eq. (8) 

Let ..., 0, Ip, 0, 0,..., 0 0 0 be a prime period two solution of the difference equa¬ 
tion Wn+i = , ^" 7 ^ -. Then 0 = , , f , , and ip = , , f, —r. This two equations lead 

«-|-i p+gw„+w„-i ^ l+pxp+q(p ^ l+pip+qi) ^ 

to the set of solutions (prime period two) except the equilibriums as {0 —)■ 0, 0 —)■ 1 — p}. 


13 



Let ..., 0, 0, 0, 0,..0 7 ^ 0 be a prime period two solution of the equation (8). We 
set 

'^n—1 
Vn = Wn 

Then the equivalent form of the difference equation (8) is 


^n+l W 
Un 


Vn+1 — 

P + qVn + Un 

Let T be the map on C x C to itself dehned by 


T 


u 


V 

U 


p-\-qv-\-u 


Then 


0 


is a hxed point of T^, the second iterate of T. 




u 


p-\-qv-\-u 


t2 


u 


p+qv+u 

\ ) 

g{u,v) 
h{u, v) 


where g{u, v) = 


p-j-qv-j-u 


and h{u, v) = 


_ p+qv-\-u 


P+9 


p-\-qv + u 


+v ■ 


Clearly the two cycle is locally asymp¬ 


totically stable when the eigenvalues of the Jacobian matrix evaluated at 
he inside the unit disk. 


0 


We have, 


where ^(0,0) 


t7'7’2 


0 


“d ‘mi-) 


(p+tp+qip) 


£( 0 , 0 ) £( 0 , 0 ) 

£( 0 , 0 ) £( 0 , 0 ) 

._ <]<t> 

{p+(t>+qp)^ 


^((h = (i+b)(p+gh) j ;v ^ _ <t>(j>+q+<l>+2qp) 

(p+0-|-ij(ji+(p+q+0)i/,-|-(jy2)2 <5t, VV^, V/ (p+0+(j0-|-(p+ij+</,)'i/,+gV’^ 


Now, set 

= + —(S = P + g 0 _ 0 (p + g + 0 + 2qi)) _ 

^ (5m ’ (5m ’ (p + 0-fg0)2 (p q. 0 + + (p + g + 0)^ g^2)2 
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_0(P + #)_ 

{p + (j) + qip) {p + (j) + q(j) + {p + q + + qij‘^Y 


A = 


For the prime period 2 solution, (j) ^ 0,'ip ^ 1 — p, x = and |A| = 0. Therefore, 


by the Linear Stability theorem (|x 
asymptotically stable if and only if 

words, the condition reduces to |p| < 1 and 
the real set up. 


< 1 + |A| < 2) the prime period 2 solution is locally 
1 


^ In other 
g| > 1 which is same condition as it was for 


< 1. It turns out that |g| > 


6 Chaotic Solutions 

This is something which is absolutely new feature of the dynamics of the difference 
equation (1) which did not arise in the real set up of the same difference equation. Com¬ 
putationally we have encountered some chaotic solution of the difference equation (8) 
for some parameter values which are given in the following Table. 1. 

The method of Lyapunov characteristic exponents serves as a useful tool to quantify 
chaos. Specihcally Lyapunav exponents measure the rates of convergence or divergence 
of nearby trajectories. Negative Lyapunov exponents indicate convergence, while pos¬ 
itive Lyapunov exponents demonstrate divergence and chaos. The magnitude of the 
Lyapunov exponent is an indicator of the time scale on which chaotic behavior can be 
predicted or transients decay for the positive and negative exponent cases respectively. 
In this present study, the largest Lyapunov exponent is calculated for a given solution 
of hnite length numerically [TO] . 

From computational evidence, it is arguable that for complex parameters p and q which 
are stated in the following table the solutions are chaotic for every initial values. 


Parameters p, q 

Interval of Lyapunav exponent 

p = (0.2037,0.5444), 
q = (0.8749,0.1210) 

(0.3215,1.6235) 

p = (0.4933,0.7018), 
q = (0.8878,0.0551) 

(1.062,2.021) 

p = (0.7840,0.4867), 
q = (0.4648,0.1313) 

(0.6256,1.314) 

p = (0.2308,0.6580), 
q = (0.5629,0.2818) 

(1.373,2.325) 


Table 1: Chaotic solutions of the equation (8) for different choice of parameters and 
initial values. 
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The chaotic trajectory plots including corresponding complex plots are given the follow¬ 
ing Fig. 1. 




Real part of The Trajectory 



Real part of The Trajectory 





Figure 1: Chaotic Trajectories of the equation (8) of four different cases as stated in 
Table 1. 


In the Fig. 1, for each of the four cases ten different initial values are taken and 
plotted in the left and in the right corresponding complex plots are given. From the 
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Fig. 1, it is evident that for the four different cases the basin of the chaotic attractor is 
neighbourhood of the centre (0, 0) of complex plane. 


7 Some Interesting Nontrivial Problems 

Open Problem 7.1. Does the dijferenee equation have higher order periodic cycle? If 
so, what is the highest periodic cycle? 

Open Problem 7.2. Find out the set of all parameters p and q for which the difference 
equation (8) has chaotic solutions. 

Open Problem 7.3. Find out the subset of the D of all possible initial values Zq and Zi 
for which the solutions of the difference equation are chaotic for any complex parameters 
p and q. Does the neighbourhood of (0, 0) is global chaotic attractor? If not, are there 
any other chaotic attractors? 


8 Future Endeavours 

In continuation of the present work the study of the difference equation Zn+i = 
where Pn, 7n, Bn and C„ are all convergent sequence of complex numbers and 
converges to a, (3, 7 , A, B and C respectively is indeed would be very interesting and 
that we would like to pursue further. Also the most generalization of the present rational 
difference equation is 

_ a + flZn-l + 'IZn-k 

A+BZn-l + CZn-k 

where I and k are delay terms and it demands similar analysis which we plan to pursue 
in near future. 
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